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Abstract 

From e + e~ collision data acquired with the CLEO detector at CESR, we observe the non-DD 
decay ^(3770) — > 7Xci with a statistical significance of 6.6 standard deviations, using the two- 
photon cascades to J/ip and J/ip -»■ £+£'. We determine a(e + e~ -> V(3770)) x #(^(3770) -> 
7Xcl ) = (18.0 ± 3.3 ± 2.5) pb and branching fraction fi(V>(3770) -► 7x c i) = (2.8 ± 0.5 ± 0.4) x 10~ 3 . 
We set 90% C.L. upper limits for the transition to Xc2 (Xco) : x B < 5.7 pb (< 282 pb) and 
B < 0.9 x 10~ 3 (< 44 x 10~ 3 ). We also determine T(^(3770) -► 7x c i)/T(V>(3770) -► ir+ir-J/ij>) = 
1.5 ± 0.3 ± 0.3 (> 1.0 at 90% C.L.), which bears upon the interpretation of X(3872). 
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Transitions from 0(3770) to other charmonium states are interesting because they test 
models of 2 3 5i — l 3 Di mixing and probe amplitudes for direct transitions from ID to 15 or 
IP states. The latter have been of considerable interest since the discovery of the narrow 
X(3872) state in tt + it~ transitions to J /if) 0, and its possible interpretation as a 1 3 D 2 
state, competing with the DD* molecule hypothesis. Measurement of hadronic transitions 
between -0(3770) and J /if) is a subject of a separate paper In this Letter, we present 
an analysis of photon transitions between -0(3770) and XcjO-P) states, followed by another 
photon transition to J/ip, with J/ip decaying to e + e _ or 

The data were acquired at a center-of-mass energy of 3773 MeV with the CLEO-c detector 
|H operating at the Cornell Electron Storage Ring (CESR), and correspond to an integrated 
luminosity of 281 pb _1 . The CLEO-c detector features a solid angle coverage of 93% for 
charged and neutral particles. The cesium iodide (Csl) calorimeter attains photon energy 
resolutions of 2.2% at E 1 = 1 GeV and 5% at 100 MeV. For the data presented here, the 
charged particle tracking system operates in a 1.0 T magnetic field along the beam axis and 
achieves a momentum resolution of 0.6% at p = 1 GeV/c. 

We select events with exactly two photons and two oppositely charged leptons. The 
leptons must have momenta of at least 1.4 GeV. We distinguish between electrons and muons 
by their energy deposition in the calorimeter. Electrons must have a high ratio of energy 
observed in the calorimeter to the momentum measured in the tracking system (E/p > 0.7). 
Muons are identified as minimum ionizing particles, thus required to leave 150 — 550 MeV 
of energy in the calorimeter. Stricter lepton identification does not reduce background in 
the final sample, since all significant background sources contain leptons. Each photon must 
have at least 60 MeV of energy and must be detected in the barrel part of the calorimeter, 
where the energy resolution is best. The invariant mass of the two photons must be at least 
3 standard deviations away from the nominal tt° or 7] mass. The total momentum of all 
photons and leptons in each event must be balanced to within 50 MeV. The invariant mass 
of the two leptons must be consistent with the J/ip mass within ±40 MeV. The measured 
recoil mass against two photons is required to be within —4 and +3 standard deviations 
from the J/ip mass. An average resolution of the recoil mass is 16 MeV. To reduce Bhabha 
background in the dielectron sample we require an average of the cosines of the angle between 
the electron direction and the direction of the electron beam and of the angle between the 
positron direction and the direction of the positron beam to be less than 0.5. The event 
selection efficiencies for -0(3770) — > iXcj, Xcj ~^ iJ/^i J/i> — > A* 4 'A* - {J/i* ~^ e + e~) events 
are 23%, 29% and 25% (13%, 17% and 15%) for the Xc2, Xa and Xco states, respectively. 

After all selection cuts, we employ kinematic fitting of events to improve resolution on the 
photon energy. We constrain the total energy and cartesian components of total momentum 
to the expected center-of-mass four-vector components, which take into account a small 
beam crossing angle. We also impose a J/ip mass constraint. No cut on confidence level of 
the kinematic fit is used, since the explicit selection cuts on the constrained quantities have 
been already employed, as described above, and because the calorimeter energy response 
function is not Gaussian. These constraints improve energy resolution for the first transition 
photon by 20%. The effect of kinematic fitting is illustrated on the CLEO-c ip(2S) data 
(1.5 x 10 6 resonant decays) in Fig. [TJ These data have clean 0(25) — > 1Xc2,i signals in 
events, which we selected with the same criteria as described above. The separation between 
these two photon lines improves after the kinematic constraints and the detector response 
function becomes Gaussian. To verify our selections and procedures, branching fractions for 
0(25) — > IXcj —> llJ/ip decays are determined from a fit to the kinematically-constrained 
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photon energy distribution (Fig. [TJ)). The normalizations, widths and positions of two 
Gaussian shapes representing large Xc2,i signals, the normalization of small Xco signal (with 
its shape fixed to the shape of the Monte Carlo distribution), and polynomial-background 
parameters float in this fit. This cross-check gives results that are within (1—2)% (relative) of 
the recently published |H analysis using different selections and signal extraction techniques. 
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FIG. 1: Energy of the lower energy photon for ip(2S) — > 7XcJ — ► JJ-J/i^, J/ip —> events in 

the CLEO-c data; (a) before and (b) after kinematic constrains on the events (see text). The solid 
line in the bottom plot represents the fit of the Xcj signals on top of barely visible polynomial 
background (dashed line). 



The photon energy distribution for the lower energy photon in the event is plotted for 
?/>(3770) — ► 7Xcj - * llJ/^i J/i> Monte Carlo data in Fig. |21 Transitions via the 
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Xc2 and Xci states produce Gaussian distributions peaked at the photon energies generated 
in -0(3770) — > 7Xc2,i decays. Transitions via the Xco state produce a broad distribution since 
the lower energy photon is usually due to the Doppler broadened Xco l-J/ip photon line, 
and sometimes due to -0(3770) — > ■yxco decay, as these two photon lines overlap each other. 

We fit the distribution observed in the data with these three signal contributions on 
top of a smooth background represented by a quadratic polynomial. The Xc2,i signals are 
represented by Gaussian peaks. The widths of the signal peaks are fixed to the values 
predicted by the Monte Carlo simulations (<te 7 = 5.1 MeV). Amplitudes of both Gaussians 
and the energy of the xa peak are free parameters in the fit. The energy of the Xc2 peak 
is constrained to be the latter minus the mass difference between these two states. The Xco 
signal shape is fixed to the Monte Carlo distribution (Fig. |2J. 

In addition to e + e~ — > ■0(3770), ?/>(3770) — > •yXcj, also e + e" — > jip(2S), 0(25') — > 7Xcj 
can contribute to the observed peaks. The cross section for the latter process peaks for 
small energies of the initial state radiation photon. Hence the produced ip(2S) mass from 
the high-mass tail of this resonance peaks at the center-of-mass energy. This makes the 
ip(2S) background indistinguishable from the ■0(3770) signal. We estimate the size of this 
background from the theoretical formulae, which fold in radiative flux, W(s,x), the Breit- 
Wigner shape of ip(2S), BW(s'), the branching ratio, Bx, for ip(2S) — > ^Xcj — ► llJ/^P 
77^+£- @ at the 0(2S) peak, and a phase-space factor, Fx(s'), rescaling the latter to the 
actually produced mass of if>(2S) at its resonance tail. Here, s is the center-of-mass energy 
(3773 MeV) squared, s' is the mass-squared with which the if)(2S) resonance is produced, 
and x is the scaled radiated energy in e + e~ — > jil)(2S), x = 1 — s'/s. Above, we have used 
the notation from Ref. 0, where the formula for W(s,x) is given and discussed in detail. 
Our selection cuts limit this radiated energy to less than 50 MeV (x < 0.027), therefore, the 
ip(2S) contribution is limited to its component which peaks near x ~ 0, where the energy 
resolution smears it to look like the ^(3770) signal. The phase-space factor Fx(s') is equal to 
(E^(s') / E^ eak ) 3 p, where Ej(s') and E^ eak are the energies of the photon in the ip(2S) — > 

7Xcj transition at the ip(2S) resonance tail (v^s 7 ~ 3773 MeV) and peak (v^s 7 = Mr), 
respectively. The ip(2S) resonance mass (Mr) and total width (Tr) in the Breit-Wigner 
formula, BW(s') = 127rr J? r ee /[(s / - M R ) 2 + M R T R ), are fixed to the world average values 
0|, while the T ee is fixed to the value recently determined by CLEO Integrating the 
theoretical cross section in the x < 0.027 range, and multiplying it by the event selection 
efficiencies given previously, we estimate that the if>(2S) background contributes 12.2, 21.1 
and 0.7 events to the Xc2, Xa and Xco peaks, respectively. The systematic uncertainty 
in these estimates is 25%. We represent these background peaks in the fit to the energy 
spectrum by the same shapes as described previously for the signal contributions with the 
amplitudes fixed to the estimated number of background events. 

The smooth background under the peaks is significantly higher in the 77e + e~ sample 
than in the 77/z + /x~ sample due to a high cross section for radiative Bhabha scattering. 
Therefore, instead of adding the photon energy distributions for these two samples, we fit 
them simultaneously, as illustrated in Fig. El The ratios of the peak amplitudes between 
the dimuon and dielectron samples are fixed to the ratios of the selection efficiencies. The 
signal shapes are constrained to be the same. The background-polynomial parameters are 
independent. 

The fitted signal amplitudes (quoted for the sum of the dimuon and dielectron samples) 
are O.Oig!o5 53 ± 10 and 22 ± 9 events for x c2 , Xci and XcO; respectively. To estimate a 
probability that the data contain no signal contribution, we also perform fits with the signal 
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FIG. 2: Energy of the lower energy photon for the simulated ^(3770) — > 7% c j — > jjJ/ip, J/ip — > 
£ + £~ events, for J = 2, 1 (solid-line histograms) and (dashed-line histogram). The vertical axis 
gives the number of detected Monte Carlo events per bin divided by the total number of generated 
events and then multiplied by a hundred. Thus, the area under each peak gives the detection 
efficiency in percent. The upper range of the horizontal axis reaches the kinematic limit. 



amplitude fixed at zero. The ratio of the fit likelihoods is transformed into the number 
of standard deviations (a) at which the null hypothesis can be excluded, which, for our 
■^(3770) — > 7Xci signal, is 6.6cr. The fitted peak energy, 253.5 ±1.2 MeV (statistical error 
only), is in excellent agreement with the 253.6 MeV value expected from the center-of-mass 
energy and the Xci mass. The data in the Xci signal region exhibit the expected peaking 
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TABLE I: Various quantities for ^(3770) — > j Xc j transitions. Efficiencies given here are averaged 
over the 77/x + \x~ and 77e + e~ channels. The upper limits are at 90% C.L. 





J = 2 


J = 1 


J = 


signal events 


o n+^- y 
u - u -o.o 


53 ± 10 


22 ±9 


efficiency (%) 


18 


23 


20 


a(e+e~ -» V(3770)) 








xfi(^(3770)^7Xcj) (pb) 


< 5.7 


18.0 ±3.3 ±2.5 


< 282 


£(</>(3770) - 1X cj) (%) 


< 0.09 


0.28 ± 0.05 ± 0.04 


< 4.4 


r(^(3770) -> 7Xc j) (keV) 


< 21 


67 ± 12 ± 12 


< 1050 



of the dilepton mass and of the two-photon recoil mass at the nominal J/i[) mass as shown 
in Fig. |U Since the statistical significances of the Xc2 and Xco contributions are 0.0 and 1.7 
standard deviations, respectively, there is no evidence for photon transitions via these states 
and we set upper limits on their rates. 

The integrated luminosity of the datasets was measured using e + e~, and 77 events 

8]; event counts were normalized with a Monte Carlo simulation based on the Babayaga 
9| event generator. The resulting systematic error in luminosity measurement is 1%. The 
systematic error in efficiency simulation is 4%. Variations in the fit range, order of the 
background polynomial, bin size and the signal width result in a variation of the Xd signal 
yield by 6%, while the systematic uncertainty in the subtraction of the ip(2S) background 
contributes 7%. An additional systematic uncertainty of 6% comes from the Xci — > iJ/ip and 
J/ip —> £ + £~ branching ratios p used in unfolding the measured rate for the ^(3770) — > 7Xci 
component. The systematic errors on the Xc2 and Xco rates are obtained in a similar way. 
To obtain upper limits, we combine statistical and systematic errors in quadrature. The 
results for a(e + e" -> ^(3770)) x £(>(3770) -> 7 x cJ ) are (18.0 ±3.3 ±2.5) pb for Xd, < 5.7 
pb (at 90% C.L.) for X c2_, and < 282 pb (at 90% C.L.) for X co- 

Using a(e + e~ — > DD) |10J for a(e + e~ — > ^(3770)), given that all measured non-DD 
decays of ^(3770) g llll | have very small cross sections, we obtain the following branching 
ratio results: 25(^(3770) -> jxa) = (2.8 ±0.5 ±0.4) x 10~ 3 , 25(^(3770) -> 7x c2 ) < 0.9 x 10~ 3 
(90% C.L.) and 23(^(3770) -> 7x c o) < 44 x lO" 3 (90% C.L.). 

We turn the branching ratio results into transition widths using r tot ('?/>(3770)) = (23.6 ± 
2.7) MeV 0. This leads to: T(^(3770) -> 7% cJ ) = (67 ± 12 ± 12) keV for X d, < 21 keV 
(90% C.L.) for X c2, and < 1.0 MeV (90% C.L.) for X co (see TableJTlfor the summary). These 
results agree well with most of the theoretical predictions [l2[ Il3l Il4j| as shown in Table El 

Combining this measurement with our determination of the 7i + 7i~J/ip rate g] we obtain 
r(^(3770) -> 7Xci)/r(^(3770) -> tt+tt- J/ip) = 1.49 ± 0.31 ± 0.26 (> 1.0 at 90% C.L.). 
The transition widths measured for ^(3770), which is predominantly the 1 3 D\ state, are 
theoretically related to the expected widths for the 1 3 -D2 state. The ratio above is expected 
to be a factor 2-3.5 larger for the 1 3 -D2 state with a mass of 3872 MeV than for the ^(3770) 
HI H Gi}. In view of the upper limit from Belle, T(X(3872) -> 7 x c i)/r(X(3872) -> 
Tr+TT-J/ip) < 0.9 (90% C.L.) Q, the 1 3 D 2 interpretation of X(3872) is strongly disfavored, 
which is also supported by other recent Belle results 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
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TABLE II: Our measurements of the photon transitions widths (statistical and systematic errors 
have been added in quadrature) compared to theoretical predictions. 





r(V(3770) -> TXcj) (keV) 




J = 2 


J= 1 


J = 


CLEO data 


< 21 


67 ± 17 


< 1050 


Rosner [12] 


24 ±4 


73 ±9 


523 ± 12 


Eichten-Lane-Quigg [131 








naive 


3.2 


183 


254 


with coupled-channels corrections 


3.9 


59 


225 


Barnes-Godfrey-Swanson [14] 








non-relativistic potential 


4.9 


125 


403 


relativistic potential 


3.3 


77 


213 



dation and the U.S. Department of Energy. 
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FIG. 3: Energy of the lower energy photon for the selected e + e~ — > 77 J/^, J/tp — > (top) 
and J/^> — > e + e~ (bottom) events at the V(3770) resonance. The solid line shows the fit. The 
dotted line shows the smooth background. The dashed line shows the total background including 
the expected background-peaks from radiatively produced tail of the ip(2S) resonance (see text). 
The latter saturates the Xc2 contribution. The excess in the Xci peak above the 2S contribution 
(dashed line) represents evidence for ^(3770) — > 7Xci transitions. 
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FIG. 4: Distributions of the J/ip mass reconstructed either as dilepton mass (left plots) or diphoton 
recoil mass (right plots) for events with the lower photon energy within ±2<r of the Xci peak. The 
cuts on both plotted quantities have been loosened to ±100 MeV to avoid selection bias on the 
displayed distributions. The points with error bars represent the data. The dashed histograms 
represent the expected amount of vp(2S) — > 7Xci background. The solid histograms represent this 
background contribution plus the ^(3770) — > 7Xci signal contribution, as simulated with Monte 
Carlo, normalized to the number of signal events determined by the fit to the photon energy 
distribution. The r y^ye + e~ data (bottom plots) have a higher level of other backgrounds and lower 
signal efficiency than the 77/i + /i~ data (top plots). 
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